Ground-state properties of few dipolar bosons in a quasi-one-dimensional harmonic trap 
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We study the ground state of few bosons with repulsive dipole-dipole interaction in a quasi-one-dimensional 
harmonic trap by means of the exact diagonalization method. Up to three interaction regimes are found depend- 
ing on the strength of the dipolar interaction and the ratio of transverse to axial oscillator lengths: a regime where 
the dipolar Bose gas resembles a system of weakly (^-interacting bosons, a second regime where the bosons are 
fermionized, and a third regime where the bosons form a Wigner crystal. In the first two regimes, the dipole- 
dipole potential can be replaced by a 5 potential. In the crystalline state, the overlap between the localized wave 
packets is strongly reduced and all the properties of the boson system equal those of its fermionic counterpart. 
The transition from the Tonks-Girardeau gas to the solidlike state is accompanied by a rapid increase of the in- 
teraction energy and a considerable change of the momentum distribution, which we trace back to the different 
short-range correlations in the two interaction regimes. 
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I. INTRODUCTION 

Ultracold atoms and molecules with large permanent dipole 
moments are currently attracting much interest, since they al- 
low realization of quantum gas systems with long-range in- 
teractions. Major steps into this direction have already been 
done. References [1, 2] reported strong dipolar effects in 
a Bose-Einstein condensate (BEC) of ^^Cr, which possesses 
a large permanent magnetic dipole moment of comparable 
strength as the usual S interaction. Even more promising 
are ultracold molecules of two different atomic species, since 
they have much larger permanent electric dipole moments 
[3]. They have already been produced by means of radio- 
frequency (rf) spectroscopy [4] and brought into the lowest 
internal vibrational ground state by means of a stimulated Ra- 
man adiabatic passage [5, 6]. What remains is to cool these 
gases down into the quantum degenerate regime. 

Many new effects have been predicted for ultracold quan- 
tum gases with dipole-dipole interactions (DDIs), which are 
based on the long range and anisotropy of the DDI. Among 
others, we mention the stabilization of a dipolar BEC in a 
pancake-shape trap [7], the roton-maxon character of the ex- 
citation spectrum [8], new exotic quantum phases in optical 
lattices [9-11], and the transfer of spin into angular momen- 
tum similar to the Einstein-de Haas effect in ferromagnets 
[12, 13]. Apart from that, the shape and strength of the inter- 
molecular interactions may be controlled by means of static 
electric and microwave fields [10], and even three-body in- 
teractions may be realized in optical lattices [11]. Moreover, 
novel quantum computation schemes are proposed with ultra- 
cold dipolar molecules [14—16]. 

In this article, we discuss different interaction regimes of 
a quasi-one-dimensional (ID) dipolar Bose gas, in which the 
permanent dipole moments are aligned by an external field, 
such that the effective ID DDI is repulsive. Thus far, theoret- 
ical studies of ID dipolar bosons found correlations [17, 18] 
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beyond those in a Tonks-Girardeau (TG) gas [19-21] and ex- 
citations [22-24], which can be described within the Luttinger 
liquid framework [25]. We show that the l/|a;p potential used 
in these references, which is valid in traps of zero transverse 
width, gives rise to an infinitely strong S peak at equal particle 
positions, even for the infinitesimal strength of the DDI, due 
to its singular behavior at x = 0. Hence, the ID Bose system 
forms a TG gas for very weak interactions [17, 18,22]. We use 
a ID DDI potential that accounts for the finite width of the trap 
and is hence finite at a; = [27]. We show that this potential 
acts like a S peak of finite strength, when the trap anisotropy is 
large. This allows for a regime of weakly (5-interacting bosons 
below the TG regime. In the different interaction regimes, we 
analyze the contributions to the total energy, the particle den- 
sity, the momentum, and occupation number distribution of 
the particles in the harmonic trap. Different features of the 
momentum distribution and the associated kinetic energy are 
related to the correlations between the particles. 



II. MODEL HAMILTONIAN 

We consider electric or magnetic dipoles, which are aligned 
in the xz-plane by an external field (see Fig. 1). The DDI 
between two point-like dipoles is modeled by 

T4d(rO = ^(l-3cos2 0,<j), (1) 

where (P is the strength of the DDI and cos 6rd = r ■ d/ [rd). 
The coupling strength between two permanent electric dipole 
moments is given by d^ = (D debye)^/(47reo), where D 
is the dipole moment in debyes and eo is the electric con- 
stant. In the case of permanent magnetic dipole moments, 
d^ — /xo.g|/i^/(47r), where /xg is the magnetic constant, 
is the Lande factor, and fis is the Bohr magneton. In the fol- 
lowing, we neglect the short-range van der Waals interaction 
[26, 27]. 

The dipoles are enclosed in a cigar-shaped harmonic trap. 
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FIG. 1 : (Color online) Dipoles d, which are oriented within the x-z 
plane and which enclose an angle Q with the x axis. 



which may be generated by a deep optical lattice [20]: 



Here, uj and wj^ are the trap frequencies of the axial and per- 
pendicular directions, respectively, and w ^ wj^. Under the 
condition that the energy per particle of the axial direction is 
much smaller than the transverse level spacing huj^, one can 
assume that the particles stay in the ground state of the trans- 
verse harmonic oscillator, exp[— (j/^ + z'^)/(2Z^)] /(Zj^Y^) 

with Zj^ ~ ^ h/{mu}±). Integration over the transverse direc- 
tions yields the effective ID DDI [27] 



with 



Kid(x) = [/ddKid(k|Ai), 
+ 3 008(26*)] 



C/dd 



(2) 



(3) 



and 



Vaa{u) = -2u + (1 + u^) e"'/2 erfc(u/\/2), (4) 




FIG. 2: (Color online) Dimensionless DDI potential Kid(|2;|//x) for 
l±_ = 0.4 1 (blue) and l± — I (red). Vdd is finite at the origin and be- 
comes more peaked for smaller Zx . At large distances, Vdd oc 



where erfc is the complementary error function. An explicit 
calculation of Eqs. (2)-(4) is done in Appendix A. The di- 
mensionless DDI potential Vdd(|a;|/Zj_) is plotted in Fig. 2. 
The second-quantized many-particle Hamiltonian is then 



aiQk , (5) 



ijkl 



where a| (oi) are bosonic creation (annihilation) operators for 
one particle in energy eigenstate (pi {x) of the axial harmonic 
oscillator and where 



dxdx' (j>^{x)(l)j{x')Vdd{\x - x'\/lj_)(t>kix)(l)i{x') 



are dimensionless interaction integrals. These integrals are in- 
tegrated numerically for different ^j^. The Hamiltonian matrix 
(5) is diagonalized in the subspace of the energetically lowest 
eigenstates of the noninteracting many-particle problem. 



III. DISCUSSION OF THE EFFECTIVE ID DDI 

We solve numerically the many-body Schrodinger equation 
with the interaction of Eqs. (2)-(4), but to gain some intuitive 
understanding let us first analyze the interaction potential. The 
DDI strength C/dd, as defined in (3), is negative for angles be- 
tween ^ 6 < f?crit. with 6'ci-it. = arccos(l/\/3) and positive 
for e'crit. < < 7r/2. For = 0, the ID DDI is maximally 
attractive with C/dd = ^'^^/(2 ^j^), and for 9 = tt/2, it is max- 
imally repulsive with C/dd ~ In the following, we 
restrict the discussion to the repulsive case C/dd > 0. 

From a Taylor expansion of (4) around infinity, one finds 
that Vdd(u) — >■ for u — > oo. Thus, for large distances 

I a; I ^ l±, the effective ID DDI is given by 

Vaaix) « U,,/\x/l\^ 

with the long-range interaction strength Uh- = 4X^Udd, where 
A — l±/l. On the other hand. Fig. 2 suggests that Vdd(a;) 
becomes a S peak for small l± . Indeed one finds 



dx-^Vdd{\x\/h) 



1, 



which leads us to the definition = Vdd(|2:|/Z_L)/(4 Zj^). 
We define the width of the function according to 



wij^ = 2v2ttIj_ w 5/j_, 
which is justified by the observation that 



dx- 



/2lT l± 



1 

'4t; 



■VM{\x\lli_) « 90%. 



Hence, the series of functions converges toward a 5 peak 
when Zj^ approaches zero; that is, 5 for Zj^ — > 0. We 

conclude that at short distances |a:;| < 2.5 Zj^ 

Vdd{x) ^ Al^UdiSix) = UJ{x), 
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with Usr = 4At/dd and 5 = S/l. It follows that the strengths of 
the short- and long-range parts of the ID DDI scale differently 
with the trap anisotropy 1/ A, 

which means that the short-range part of the interaction is 
strongly enhanced by a tight transverse confinement compared 
to the long-range part. Particularly in the limit A = 0, the ef- 
fective ID DDI is not only given by Uir/\x/l\'^ at a; 7^ 0, but 
there is an additional infinitely strong 5 peak at x = 0. Con- 
sider, for example, the situation of a rather small Uji- = 0.1 fiuj 
and a trap anisotropy of 1/A = 10. Then, the strength of the 
S peak is C/s, = 10 hu, which is already so large that bosons 
form a TG gas. 

This decomposition of the ID DDI is valid as long as the 
ground state is a slowly varying function along the axial di- 
rection on the transverse length scale l± . In the following we 
fix the trap anisotropy according to 1/A — 50 so that this 
condition is always fulfilled. This corresponds to the max- 
imal trap anisotropy, which has been achieved in the exper- 
iment of Kinoshita et al. [20], where uj = 27r27.5Hz and 
max(wj^) = 27r 70.7kHz leading to max(l/A) = 50.7. 

In Sec. IV, we perform a sweep of U^, which can be done 
by changing the angle 6 between arccos(l/'\/3) and tt/2. 
As mentioned before, we neglected the short-range van der 
Waals interaction, since we are interested in a study of the 
effect of the ID DDI alone. These relations, which we de- 
rived within the single-mode approximation, are true as long 
as ^dd > l± [27], where ^dd = dPm/h^. By taking '^OR^'^Rb 
polar molecules with an electric dipole moment of 0.6 debye 
and cjj^ = 27r 70.7kHz, the condition is fulfilled for 6* > 1. 
However, weak dipole moments may largely influence the ID 
scattering in the regime ^ Zdd, which was shown in a the- 
ory that accounts for the short-range van der Waals interac- 
tion and goes beyond the single-mode approximation [27]. To 
study the effect of the ID DDI alone in that regime, one may 
tune giD to zero using a Feshbach resonance for fixed 6 and 
and change w. An increase of w, for example, would de- 
crease t/sr/ {huj) and increase U\^/ ihiu). 



IV. RESULTS 

Let us now turn to our results for four particles obtained by 
means of a numerical diagonalization of (5). Figure 3 shows 
the contributions to the total energy as a function of C/i, (lower 
scale) and Usi (upper scale), respectively, in a double logarith- 
mic plot for a trap anisotropy of 1/A = 50. Both scales are 
related to each other through L/ji = Uh/X^. Shown are the 
energies of bosons with DDI (solid) and 6 interaction (dash- 
dotted) and of fermions with DDI (dotted). One sees that all 
energy contributions of bosons with DDI agree with those of 
(5-interacting bosons for U\r < 0.1 hw; that is, in this region 
the bosons feel only the short-range S part of the DDI. Hence, 
it is more illuminative to use the upper scale here. The transi- 
tion behavior in this interaction regime has been discussed in 
Ref [28] (for other trap geometries, see [29, 30]): The system 
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FIG. 3: (Color online) Contributions to the total energy of four par- 
ticles as a function of [/i, (lower scale) and f/sr (upper scale), respec- 
tively, in a double logarithmic plot for a trap anisotropy of 1/A = 50. 
Dipolar bosons (solid) are compared to 5-interacting bosons (dash- 
dotted) and dipolar fermions (dotted). The shaded background marks 
the TG regime, where the bosons fermionize. 



evolves from a weakly interacting quasi-BEC via an interme- 
diate regime to a TG gas. Already at J7si- ~ 10 huj, the bosons 
are fermionized, which is indicated by a saturation of the total 
energy. In the shaded region, the system does not react to a 
further increase of the interaction strength, which shows that 
the bosons do not feel the long-range part of the DDI, which is 
of the order of at distance I. However, when the DDI is in- 
creased above the critical value C/i, > 0.1 huj, the fermionized 
bosons are further pushed apart from each other by the long- 
range l/|a;p tail of the DDI. An obvious signature of these 
beyond-TG correlations is the rapid increase of the interac- 
tion energy, indicating that the l/jxj^-tail of the DDI has sig- 
nificant overlap with the many-body wave function. Clearly, 
since the bosons are fermionized, all the energy contributions 
coincide with those of fermions with DDI in this region. 

Another interesting aspect concerns the kinetic energy, 
which decreases in the l/|a;p-tail-dominated regime, while 
all the other energy contributions grow with increasing repul- 
sion. We attribute this behavior to the short-range correlations 
of the wave function, which can be approximated by [24] 



,Xn) oc 



ni 

1<J 



,1/K 



l/m 



(6) 



In the TG regime, the Luttinger exponent is K ~ 1, which 
leads to a large gradient of the many-body wave function at 
Xi = Xj. This gives rise to a rather large kinetic energy 
in these regions of the configuration space. With increas- 
ing long-range interactions, the Luttinger exponent decreases, 
K < 1 [18], which diminishes the gradient and hence the 
kinetic energy at Xi — Xj . For even stronger long-range inter- 
actions, the many-body wave function can alternatively be ap- 
proximated by localized Gaussian-like wave packets [17, 24]. 
As can be seen in Fig. 4, the repulsion between the bosons 
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FIG. 4: (Color online) Density of four particles for different interac- 
tion strengths. With increasing repulsion, the particles localize due 
to the long range of the interaction. For the shown values of C/ir, 
the densities of fermions and bosons are equal, since the bosons are 
fermionized. 



mainly increases the distance between the center points of the 
wave packets, which has no influence on the kinetic energy. 
Hence, we expect a saturation of the kinetic energy for very 
strong repulsion, which eventually becomes negligible com- 
pared to the other energy contributions; see Fig. 3. 

The main message of Fig. 3 is the distinction of three in- 
teraction regimes. In the left region, the bosons feel a 5 po- 
tential of finite strength; in the middle region, the bosons feel 
an infinitely strong 6 potential and thus fermionize; and in 
the right region, the system is dominated by the long-range 
l/|a;|^ part of the DDI. The right boundary of the middle re- 
gion is independent of the trap anisotropy 1/A, but the left 
boundary moves to the left when 1/A is increased, and in the 
limit 1/A = oo, the left region is absent. On the other hand, 
for very low trap anisotropics, the width of the middle region 
shrinks to zero. In Fermi systems, the left interaction regime 
is absent, since fermions do not feel the 6 part of the DDI. 

Since the two (5-interaction-dominated regimes have al- 
ready been discussed in Ref. [28], we concentrate in the fol- 
lowing on the 1/ 1 a;| '^-tail-dominated regime. Figure 4 shows 
localization of the particles with increasing long-range inter- 
action. In the TG regime, the particles are rather close to- 
gether, and the minor oscillations visible in the top left den- 
sity of Fig. 4 disappear for large particle numbers N . With 
increasing long-range interaction, the density resembles four 
overlapping localized wave packets, which move apart from 
each other The equilibrium positions of the wave packets 
minimize the potential and interaction energy (the kinetic en- 
ergy is negligible). Occurrence of this quasiordered state was 
identified in the homogeneous system by means of the corre- 
lation function [17] and the static structure factor [18]. 

Particularly interesting is the momentum distribution 
(Fig. 5), which is different for bosons and fermions in the TG 
regime (in contrast to the energies and densities). One sees 
that both distributions become equal for strong long-range re- 
pulsion. This reveals that the statistics of the particles be- 
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FIG. 5: (Color online) Momentum distribution of four particles for 
different strengths of Uh- The top left figure shows the momen- 
tum distribution of noninteracting fermions and fermionized bosons, 
which are different. For very strong repulsion, the distributions of 
fermions and bosons become equal. For = 32 hu, the momentum 
distributions are the Fourier transform of the localized wave packets 
of the bottom right particle density of Fig. 4. 



come unimportant if there is no significant overlap between 
the wave packets of the individual particles. 

In the TG regime (Fig. 5, top left), the bosonic momen- 
tum distribution exhibits a high zero-momentum peak and 
long-range high-momentum tails (this is more clearly visi- 
ble in Ref. [28]). The high-momentum tails originate from 
the cusps in the TG wave function at Xi — Xj [31, 32]. 
Note that although the bosonic and fermionic distributions 
are markedly different in the TG regime, their kinetic energy 
-E'kin. — J dpp^ p{p) / {2in) is equal. This shows that the high- 
momentum tails (and hence the short-range correlations) con- 
tribute significantly to the kinetic energy. 

With increasing repulsion, the zero-momentum peak of 
the bosonic distribution significantly decreases and the high- 
momentum tails vanish, which indicates a redistribution from 
low and high toward medium momenta. Finally, for very 
strong repulsion (Fig. 5, bottom right), the bosonic and 
fermionic distributions are only marginally different. The 
broad Gaussian distribution of the crystalline state is essen- 
tially the Fourier transform of the wave packets of Fig. 4, hke 
in the Mott insulator phase in an optical lattice [33]. 

We close our discussion with Fig. 6, which shows the oc- 
cupation number distribution of the ID harmonic oscillator 
states. As in the case of the momentum distribution, one sees 
that the bosonic and fermionic distributions become equal for 
strong long-range repulsion. For Ui^ = 32 hut, the distribution 
exhibits two features: a sharp peak with a maximum popula- 
tion of the first excited harmonic oscillator state and a broad 
distribution over the states i = 3 — 14 with a maximum at 
i — 7. The broad maximum shifts to larger i when is 
further increased. 
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FIG. 6: (Color online) Occupation number distribution of four par- 
ticles in tlie ID harmonic oscillator states for different values of 
U]x. The top left figure shows clear differences between fermion- 
ized bosons and noninteracting fermions. With increasing repulsion, 
the distributions of fermions and bosons become equal. 



V. CONCLUSIONS 

We analyzed several microscopic ground-state properties of 
few dipolar bosons in a ID harmonic trap. The interaction be- 
tween the bosons was modeled by a ID potential, which was 
obtained from the three-dimensional DDI potential by an inte- 
gration over the transverse directions within the single-mode 
approximation. For large trap anisotropics, the ID DDI acts 
on the ground state like the sum of a shoit-range S and a long- 
range l/|a;p potential. Depending on the relative strength 
of both contributions, the system forms a weakly interacting 
quasi-BEC, a TG gas, or a solidlike state of particles, which 
are localized due to the strong long-range repulsion. 

The transition from the TG gas to the solidlike state was 
clearly visible in the momentum distribution. With increas- 
ing interaction, one observes a redistribution from low and 
high toward medium momenta. The disappearance of the 
high-momentum tails (and hence the short-range correlations) 
is responsible for the decrease of the kinetic energy. For 
very strong repulsion, the momentum distribution has a broad 
Gaussian shape and is essentially the Fourier transform of the 
localized wave packets. In contrast to the kinetic energy, the 
interaction and the potential energy grow rapidly with increas- 
ing repulsion in the solidlike state. 

We remark, that for weak dipolar interactions one has to 
consider scattering effects beyond the single-mode approxi- 
mation [27]; see the discussion at the end of Sec. III. Hence, 
one may use Feshbach resonances in order to enter the regime 
of weak interactions in a controlled way. 
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Appendix A: Calculation of the ID DDI 

It is assumed that the particles reside in the ground state of 
the transverse directions. The effective ID DDI is given by 



Vdd(a;i ~ X2) = J dyidy2dzidz2Vdi{ri - r2) 

X 't>l{vi)<Pl{y2)4>l{zi)cl>l{z2), 



with (j)Q{u) = exp[-u2/(2/^)]/^/j_V7r and Vdd(?^i - fa) 
given by Eq. (1). We introduce relative and center-of-mass 
coordinates fi/2 ~ R±r/2 and perform the integration over 
Y and Z: 

= ^J dydzV,,{r)e-<^y'+^'^/'^^''^\ 

In cylindrical coordinates r — {x,y, z) — (x, p cos cj), p sin (j>) 
the integral becomes 

^dd(2:) = ^ J d4>dpVM{x,pA)pe~'''^''''^\ 



where Vdd(2^, p^ 4>) is given by 

d^ 



Vii{x,p, 0) 



with 



cos Urd 



d X cos 6 + p sin </> sin 



rd 



y^X^ + p^ 



Here, we assumed that the dipoles lie in the x — z plane (see 
Fig. 1), so that d/d — (cos 9, 0, sin 9). Integration of cos^ 9rd 
over (j) gives 

r2Tr 



) cos tird 



Hence, we obtain 

1 '■2- 



{2x^ cos^e + p^ sin^ 6*). 



p2 - 2a;2 



2j^l2 I d(j3VdA{x,p,4>) ^ Ai- ^_ _5 



yjx'^ + p^ 



with Ai = [1 + 3 cos(26')]/(4Z5^). It remains to perform the 
integration over p 



We substitute = x'^ + p^ (— udu = pdp) and obtain 

POO 



\x\ 
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Next, we set v = u/{^/2l±), which leads to 

Vkd(x) = Aae" / dve — — | , 

with a = \x\/{V2li_) and A2 = (f[l + 3cos{29)]/{8V^ll^). 
From an integration by parts, we obtain the recurrence relation 

"J 2 



dv- 



(n — l)a 



which is valid for n — 2,4,6, . . . . Since /q = Y^erfc(a)/2, 
the integrals I2 and I4 are given by 



a 



/tt erfc(a) 



and 



3a3 3 



/7rerfc(Q;) 



This is inserted into 
which becomes 



/2/2 



;]^e" erfc(a) 



2V^a;2e""erfc( a] 



2x' 



Using the definition of a, one sees, that the first and the third 
term cancel each other, and one obtains 



Vm{x) = A2V2 2l_i_\x\ - V2n{ll + x'^)e°' erfc(a) 



which equals (2)-(4). 
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